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In the first part of this lecture, the X/N expansion technique is illustrated 
for the case of the large- N sigma model. In large- N gauge theories, the X/N 
expansion is tantamount to sorting the Feynman diagrams according to their 
degree of planarity, that is, the minimal genus of the plane onto which the diagram 
can be mapped without any crossings. This holds both for the usual perturbative 
expansion with respect to powers of g 2 = g 2 N, as well as for the expansion of 
lattice theories in positive powers of X/g 2 . If there were no renormalization effects, 
the g expansion would have a finite radius of convergence. 

The zero-dimensional theory can be used for counting planar diagrams. It 
can be solved explicitly, so that the generating function for the number of 
diagrams with given 3-vertices and 4-vertices, can be derived exactly. This can be 
done for various kinds of Feynman diagrams. We end with some remarks about 
planar renormalization. 



1 Introduction: The Large N Sigma Model 

A simple example of a theory where new and interesting information can be 
obtained by studying it in the limit where the number N of fundamental field 
degrees of freedom becomes large, is the Large N Sigma Model 1 . Let <j>i(x) 
be a set of N scalar fields, i — 1, ■ ■ ■ , N , and consider the Lagrangian of a 
renormalizable theory: 

N / N \ 2 

£ = -IE[(^) 2 +^]-i a > w 

i=l \i=l I 

where A is a coupling constant, usually taken to be positive. 
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The Feynman diagrams of the theory consist of propagators, 



I 





(2) 



k 2 + to 2 — ie 



to be connected by vertices: 



X 



./ 



k 



(3) 



We write this vertex as the sum of three vertices, to be indicated the following 
way: 



we get the same Feynman rules [since integrating out A yields back Eq. (1)], 
but now the double line is the A propagator, —A, and the A lines form 3- 
vertices with the <f> lines. It is now relatively easy to integrate out the field <f> 
in (5). We must remember that A is x dependent, so we do this perturbatively 
in the A field. This way, we get effective interactions in the A field, the 
interaction being described by 'vertices' formed by closed <fi lines, see Fig. 1. 
Notice that the set of Feynman rules obtained this way is easily seen to be 
the same as the ones generated by Eqs. (l)-(4), as of course it should be. 



Since (5) is only quadratic in 4> the effective vertices only consist of single <j) 
loops, and so they are easy to calculate. 




Rewriting the Lagrangian (1), using an auxiliary field A(x), as 





o 



Figure 1. Effective vertex for A field. 
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Since there are N fields <f>, each of these effective vertices go with a factor 
N. It is now instructive to replace the field A by A/yN, which amounts to a 
propagator —XN for A, whereas the i/-point effective vertices each receive an 
extra factor N~ v / 2 . Together with the factor N mentioned above, we find that 
the effective 2 point vertex goes with a factor N° (so it is iV-independent), 
the 3-vertices go with TV -1 / 2 , etc. The tadpole diagram (the effective vertex 
with v = 1) must be assumed to vanish, after application of a normal ordering 
procedure. It would be proportional to N +1 / 2 . 

If now we choose the limit N — > oo, while A = AiV is kept fixed, then 
we see that a perturbative theory emerges, where all effective vertices with 
v external lines are suppressed by factors N 1 ~ v / 2 , so that 1/s/N acts as a 
coupling constant. Only the 2-point vertex is non-perturbative, so it has 
to be computed exactly, and added to the 'bare' term —A to obtain the A 
propagator. Fortunately, the 2-point vertex correction F(m,q), where q is 
the momentum of the A field and m is the <f) mass, is easy to compute. In n 
spacetime dimensions, 

= C(n) / dx[m 2 + q 2 x(l-x)]2~ 2 , (6) 
J o 

where C{n) = \{Ait)- n / 2 T(2 - §). 

This has to be added to the bare A propagator term, so the series of 
diagrams contributing to the dressed propagator amounts to: 

= + o + 

This is the exact propagator to lowest order in 1/N. The effective Feynman 
diagrams describing the 1/N corrections are sketched in Fig. 2. In many 
respects, the 1/N expansion is like the coupling constant expansion in a field 
theory Thus, the effective propagators in diagrams such as Fig. 2 can form 
closed loops. 

At n = 4, the function C(n) has a pole, which has to be absorbed into 
a redefinition of f/A, as usual. At n — 6 or higher, there are divergences 
which require more troublesome subtractions, since there the theory is non- 
renormalizable. It is instructive to look at the q dependence of this effective 
propagator. At n < 4, this propagator is regular at high ±g 2 , since F tends to 
zero. F is complex if q 2 < —Am 2 . If m is low enough, and/or A large enough, 
then a new pole may emerge: an (unstable) bound state of two 4> particles. 



thooft: submitted to World Scientific on February 1, 2008 3 




Figure 2. Effective Feynman graphs for the N — * oo Sigma Model. 



If, however, n > 4, the function F diverges at large ±q 2 . Note that the sign 
of F after the required subtraction flips, so that for large q 2 F behaves as 

F(m, g 2 )^-3^1ogg 2 . (8) 

Consequently, in Eq. (7) a pole at some large positive value of q 2 is inevitable. 
Thus, we have a tachyonic pole in the propagator: the theory has a Lan- 
dau ghost. This ghost is directly related to the fact that the theory is not 
asymptotically free. At very large momenta, the effective quartic coupling is 
negative, which implies a fundamental instability of the theory. This is an 
important lesson from the study of the l/N expansion of the Sigma Model: 
lack of asymptotic freedom in the ultraviolet makes the theory inconsistent. 
There is no tendency to go to an ultraviolet fixed point. The theory does 
not show any inclination to "cure itself." In the literature, we often see the 
belief that models with a decent-looking Lagrangian are well-defined, "if only 
we knew how to do the mathematics." Here, we happen to be able to do the 
mathematics, and find the model to be inconsistent at n > 4 due to its bad 
ultraviolet behaviour, just like any other non-renormalizable system. 
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A propagator ghost propagator quar k lines 



vertices 



Figure 3. Fcynman rules for U(N) gauge theories. 



2 1/N for Gauge Theories 

Many gauge theories are asymptotically free 2 , so here a large N expansion 
could be much more promising. The large N expansion for gauge theories can 
be formulated for all large Lie groups, and the final results are all very similar. 
Let us first consider the gauge group U(N). The vector fields are Hermitean 
matrices A1-, where both indices i and j run from 1 to N. Writing 



F„ 



F) A l 



d v Ay+ig [Aft, A v \\ 



(9) 



we have 



= -Tr(j(^A„) 2 - ^A^A^+igd^A^, A v ] - \g 2 [A„, A v \ 

-ipii^dfj, + igj fl A li +m^,)ip. (10) 

A transparent notation for the Feynman diagrams 3 is one where the propa- 
gators are represented by double lines, each line indicating how the index is 
transported, see Fig. 3. The mathematical details of these rules depend on 
the gauge fixing procedure, which we shall not further expand upon. 

Everywhere in a Feynman diagram where an index loop closes, we get 
a factor N from counting the values that the index can have. This is the 
only way the dimension of the gauge group enters. Now consider a general 
diagram with V3 3-point vertices, V4 4-point vertices, P propagators and F 
index loops. Assume, temporarily, no external lines (a gauge-invariant cou- 
pling to an external source can for instance be introduced by means of a 2 
quark vertex insertion, Jipvp; this then is a 2 point vertex, counted by V2). 
Write V = V 2 +V 3 + V A . 
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First, by drawing dots at both ends of every propagator, we see immedi- 
ately by counting dots that 

2P = Y j nV n . (11) 

n 

We now turn the diagram into a polyhedron by defining every index loop to 
be a facet. The facets form a polyhedron if there are no quark loops. Then, 
we have Euler's theorem for polyhedra: 

F-P + V = 2-2H. (12) 

Here, H is the genus, or number of 'holes' in the polyhedron. This theorem 
is easily proven by induction. Since all 3-point vertices in Fig. 3 come with a 
factor g, and all 4-point vertices with a factor g 2 , the g and N dependence of 
the amplitude is seen to be 3 

g V 3+ 2V iN F = g 2P-2V N F = {g 2 N) F+2H-2 N 2-2H ^ 

In the limit where N — * oo, while g 2 = g 2 N is kept fixed, we see that the 
N dependence of the diagrams is N 2 ~ 2H , so that diagrams without holes 
(H = 0) dominate. 

If quarks are present, we simply declare the quark loops to be facets as 
well, but they lack a factor N. Therefore, we also get a factor , where 
Q is the number of quark loops. Thus, the l/N expansion is found to be an 
expansion with respect to the number of quark loops and the number of holes 
in the polyhedra. In many respects, this expansion resembles an expansion of 
a string theory in powers of the string coupling constant. 

For other gauge groups, the result is very similar. If we take SU(N) 
instead of U(N), we are essentially subtracting a U(l) component, whose 
coupling is g = g/y/N, so even at large g, the 17(1) subtractions are still 
perturbative. For SO(N), we have an additional reality constraint on the 
fields A: 

K = -Hi - • (W) 

This means that the arrows disappear from the double line propagators in 
Fig. 2. The surfaces of the planar diagrams are now non-oriented, so that 
non-orientable diagrams are allowed: Klein bottles and Mobius strips. They, 
however, have H > or Q > 0, and hence they are also suppressed at large 
N, so that the dominating diagrams are essentially the same as the ones for 
U(N) and SU(N). 
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3 1/N for Theories on a Lattice. 



It is instructive to consider the 1/N expansion for lattice gauge theories . 
Consider a lattice with lattice sites, connected by lattice links. The gauge 
fields live on the links, the quark fields live on the sites. On a link with length 
a, connecting the sites x and x + e p (where is a unit vector of length a on 
the lattice), we write 

Un(x) d = exp(iagAf,). (15) 
The gauge-invariant Wilson action (in n dimensions) is 

n — 4 

S = \f E Tr [M*)Uu{x + eJUfrx + e v )Ut{x)] . (16) 

The two indices of the field variable U^^x) refer to gauge transformations in x 
and x + e^, respectively, so these are two different indices, to be distinguished 
by underlining one of them. Since U(x) is unitary, we have 

uM k = Sf; Ul j U 3 k = 6t, (17) 

from which we derive that the 'functional' integral over the U values will 
yield 5 : 

(UM e ) = , (18) 

[simply by observing that (18) is the only term with the required gauge in- 
variance, so all that has to be determined is the prefactor]. Similarly, we have 
only four possible terms for the quartic averaging expression. Two of these 
turn out to be leading: 

(uMu^ui) = JL (^f of + 121) + o(n- 3 ) ■ (19) 

The effect of these expressions is that when we do the functional integral, 
simply by inserting the complete expansion of the exponential of the action S, 
at every link every factor U has to be paired with a factor W . Every such pair 
produces a factor 1/N. As in the continuum case, we construct a polyhedron 
from the elementary plaquettes. Let F be the number of plaquettes (the 
number of terms of S in our 'diagram', coming from expanding e tS ), let P be 
the number of pairings, and V be the number of indices that end up being 
summed. We find a factor 1/N for each pairing, a factor N for each index, 
which amounts to a factor N for each point on our polyhedron, and a factor 
1/g 2 for every plaquette (apart from a harmless numerical coefficient). 
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All in all, the N and g dependence arises as a factor 



N v (l/N) p (g 2 )- F = (g 2 N)- F N v - p+F = (g 2 )~ F N- 



2-2H 



(20) 



where again Euler was applied, and g 2 = g 2 N. So again we find that the 
diagrams only depend on g, and the N dependence of the diagrams is as in 
the continuum theory. This result is not quite trivial, since we are looking 
at a different corner of the theory, and the transition from the small g to the 
large g limit need not commute with the large N limit, but apparently it does. 

4 Divergence of Perturbation Expansion 

Gauge theories can now be expanded in terms of several parameters: we 
have the g expansion, the 1/g expansion and the 1/N expansion. Of these, 
only 1/N acts as a genuine interaction parameter, describing the strength 
of the interactions among mesons. Physical arguments would consequently 
suggest that the 1/N expansion produces terms of the form N~ k where at 
will diverge faster than exponentially for large k values, typically containing 
k\ terms. In contrast, the g expansion generates only planar diagrams. In 
principle, one could expect better convergence for this expansion. Indeed, the 
number Rk of planar diagrams with k loops will be bounded by an expression 
of the form C k for some finite constant C. 

Unfortunately, even the g expansion will generate factors fc!, not due to 
the number of diagrams being large, but because of rcnormalization effects. 
These effects already occur in the large N sigma model of Sect. 1. Consider 
the A expansion of the propagators of Figure 2. The expansion term going as 
A fe+1 contains exactly k loops. Each of these loops generates an expression 
F(m,q), which, for large q, diverges roughly like log(g 2 + m 2 ). Imagine this 
propagator itself being part of a closed loop, and imagine that a sufficient 
number of subtractions has been carried out in order to make the integral 
convergent, so that some renormalized amplitude results. We see that the k th 
term requires the calculation of an integral of the form 




(21) 



Writing q 2 + m 2 w m 2 e x , we see that integrals of the form 




(22) 



arise. 
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Such k\ terms being very similar to the ones generated by instanton effects 
in quantum field theory, this effect was dubbed 'renormalon'. The renormalon 
is closely associated with the Landau ghost. Indeed, replacing the coupling 
constant A by a running coupling constant, we notice that the Landau ghost 
mass increases exponentially with the inverse coupling constant: 

Wgh oc e 16 ^ . (23) 

In some very special cases, the Landau ghost can be avoided altogether. 
This is if the theory is both asymptotically free in the ultraviolet and infrared 
convergent due to mass terms 6 . Such theories can be constructed, making 
judicious use of gauge vector, scalar and spinor fields. Integrals of the type 
(22) do still occur in that case, but the perturbation expansion can then be 
Borel resummcd 6 . 

QCD is far more divergent than that. Here, we have no mass term to 
curb the infrared divergences. Instanton effects are to be expected at finite 
N, and possibly even at infinite N, and infrared renormalons invalidate simple 
attempts to prove convergence of Borel resummation 7 . Simple physical con- 
siderations, however, could be used to argue that unique procedures should 
exist to resum the planar diagrams. 

5 Counting Planar Diagrams 

Doing such resummations analytically is way beyond our present capabilities. 
What we can do is carry out such summations in a simple, zero-dimensional 
theory. This result will give us an analytic method to count planar diagrams. 
The 'field theory' is described by the action 8 ' 9 

S(M) = Tr {-\M 2 + i 9 M 3 + ±AM 4 ) , (24) 

where M is an N x iV-dimensional matrix, in the limit 

N^oo, g, A^O, Ng 2 =g 2 and N\ = \ fixed. (25) 

Henceforth, the tilde (~) will be omitted. Precise calculations of the gen- 
erating functions for the number of planar diagrams, with various types of 
restrictions on them, were carried out in Ref. 10 . Different cases were stud- 
ied, all related to one another by exact equations. Part of that paper is based 
on pioneering work by Koplik, Neveu and Nussinov 8 — which in turn makes 
use of earlier work by Tutte 9 — and on the work by Brezin, Itzykson, Parisi 
and Zuber 11 . The latter apply matrix theory to do the integration with the 
action (1). It is instructive instead to work directly from the functional equa- 
tions, as these will be easier to handle in the QCD case, and they are also more 
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transparent in diagrammatic approaches. The relations are read off directly 
from the diagrams. 

A delicate problem then is the choice of boundary conditions for these 
equations. They can be derived by carefully considering the holomorphic 
structure that the generating functions are required to have. Once this is 
understood, a fundamental solution is obtained for the generating function 
describing the numbers of all planar diagrams for all multiparticle connected 
Green functions, with given numbers V3 of three-point vertices, V4 four-point 
vertices, and E external lines. 

We denote the number of all connected planar diagrams by N^e,v 3 ,Va)- 
The generating function F(z,g. A) is defined by 



F(z,g,X) 



OO, OO, OO 

E=l, 
v 3 =o, V 4 =0 



(26) 



Here, it is for technical reasons that we start counting at E = 1. 



The recursion equation for the generating function F = F{z, <?, A), is 



F= z+ g(F 



F -F< 



(0) 



+ A [F 



2F- 



F-F, 



(0) 



F - F(o) - F (i) z 



F, 



F-F, 



(0)- 



(0) 



(27) 



Here, F( ) = F(Q,g, A) is the tadpole diagram (E ~ 1), and F^ = 
-^F(z, g, X)\ z=0 is the propagator (E — 2). Written diagrammatically, we 
have 



+ 



F - 



+ 



(a) 



F = 



F = + 



F = 



(b) 



(c) 



(d) 



F = 




(e) (f) (g) (h) 

Note that inside the 1-loop terms the tadpole diagrams had to be removed 
by hand, and in the 2-loop term the propagator had to be removed. Such 
equations can be solved! There are three methods: 



thooft: submitted to World Scientific on February 1, 2008 10 



i. The first terms in these lists of expressions can be found numerically. 
We could then try to guess the generic expressions for N^e-i,v 3 ,V4); ano - 
check these expressions in the closed equation (27) or (28). This works 
only in the very simplest of cases. 

ii. One can study the holomorphic structure of Fa(z,g, A) near the origin, 
g = X = 0. Demanding the correct non-singular behaviour at the origin 
provides the unique boundary condition. 

Hi. One can return to the definition of the generating functions using the 
integral expression for large matrices, and take the limit for the infinite 
matrices M. The integrand contains only one matrix, which can be con- 
veniently diagonalized (in case of two or more fields, or in dimensions 
higher than zero, simultaneous diagonalization of the matrices M a (x) is 
impossible, which leads to considerable complications). This turns out 
to be the most elegant method. 

After some work, we found that all approaches lead to consistent results. One 
first derives the values of and Fm, after which the values for F(z,g, A) 
follow by formally solving Eq. (27). 

The expressions found for F(0,g, A) turn out to be singular on a curve 
in the g-X plane. This gives us the critical values of g and A, in terms of a 
parameter t that parametrizes this curve: 



g c = 2t 2 (3-t)(6-6t-3i 2 + 2i 3 ) 1 / 2 (12-4< 3 + i 4 )- 3 / 2 ; 

X c = (2 + 2t - t 2 ) (6 - 6t - 'it 2 + 2< 3 ) (12 - 4i 3 + i 4 )~ 2 . (29) 



All real values for g are allowed, but A must be positive. The curve defined 
by Eqs. (29) in the g-X plane then looks as in Fig. 3. Note that it is the g 
dependence, not the g 2 dependence that we are looking at. The cusp at g = 
is a genuine singularity of the curve (29). 

Next, we can make restrictions on the occurrence of insertions within the 
diagrams. This we indicate by using a subscript a, taking values between 
and 5: 



a = 



no further restrictions. 



(30) 



a = 1 



no tadpoles: 




(31) 




(32) 



a 



= 3 



also no self-energies: 




(33) 
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(0,X C ) 








8 



Figure 4. Line of critical A and g values to a high precision. The lines are not exactly 
straight, and in fact form a single curve with a cusp singularity. The corners are given by 

the values A c = j~ and g c = l/yl^VS. 



a = 4 : also no dressed 3- vertices: 



a = 5 : also no dressed 4- vertices: 



(34) 
(35) 



In each of these cases, there is a domain in the g-X plane where the 
zero-dimensional theory converges. The location of the boundaries of these 
domains then tells us how the number of diagrams of the specific types increase 
at high orders. For example, the last curve, of the case a = 5, is given by the 
parametric equations 

c 2t 2 (6-t 2 ) 2 (6-3t 2 + t 3 ) 

95 ~ (24 - 72 1 - 12 t 2 + 56 t 3 - 10 t A - 10 1 5 + 3 i 6 ) 3 / 2 ' 

A5 = (3456 - 31104 1 + 25920 1 2 + 55296 1 3 - 180576 1 4 - 44064 1 5 

+ 247824 1 6 - 27936 1 7 - 147672 1 8 + 52760 i 9 + 38076 1 10 - 24432 1 11 



1762 1 12 + 4506 1 13 - 921 1 14 - 202 1 15 + 114 1 16 - 18 1 17 + t 18 ) 



x (24 - 72t- 12^ + 56^- 10t 4 - 10r + 3t b ) 



6\-3 



(36) 



These, as well as the other expressions that were derived in Ref. , are 
all exact. The domains of convergence for the various cases are displayed in 
Fig. 5. The last boundary, that of region # 5, is the one given by Eqs. (36). 
Figures 4 and 5 show the actual domains with great precision. All curves have 
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Figure 5. Critical regions after imposing different constraints. 



cusps at g = 0. Many of our results were obtained and/or checked using the 
computer program Mathematica. 

6 The g 2 N Expansion 

We conclude with a short discussion on constructing perturbative planar QCD 
amplitudes at all orders without divergences. 

A renormalization scheme that appears to be not so well-known 6 can be 
set up in the following way. It works particularly well for planar field theories 
in 4 space-time dimensions, but it also applies to several other quantum field 
theories, notably QED. Presumably, QCD with iV c = 3 can also be covered 
along these lines, but some technical details have not been worked out to my 
knowledge. Our scheme consists of first collecting all one-particle irreducible 
2-, 3- and 4-point diagrams, and formally considering the non-local quartic 
effective action generated by these diagrams. 

Next, consider all diagrams using the Feynman rules derived from this 
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action, but with the limitation that only those diagrams that are absolutely 
ultraviolet convergent are included. No superficially ultraviolet divergent sub- 
graph is accepted. To be precise, we omit all diagrams with 4 or less external 
lines, as well as all diagrams containing any non-trivial irreducible subdiagram 
with 4 or less external lines. They are the diagrams of type 5 of the previous 
section. 

Clearly, one expects that the 1PI subgraphs with 4 or less external lines 
have already been taken care of by our use of the quartic effective action in- 
stead of the bare Lagrangian. The important issue to be addresed is, whether 
the counting of all diagrams was done correctly so as to obtain the required 
physical amplitudes. But this is not difficult to prove: 

Theorem: The above procedure correctly reproduces the complete ampli- 
tudes for the original theory. No diagrams are over-counted or under-counted. 

The proof of the theorem is to be found in Refs 6 ' 10 . One deduces that 
if all one-particle irreducible 2-, 3-, and 4-point vertices are known, all am- 
plitudes can be derived without encountering any divergent (sub-)graph(s). 

Subsequently, what has to be done to complete a perturbative computa- 
tional scheme, is to establish an algorithm to compute the irreducible 2-, 3-, 
and 4-point vertices (and, if they occur, the tadpole diagrams as well). Ac- 
tually, this is simple. Consider a 4-point 1PI diagram T(pi_,p2,p$,Pi). Here 
p^ are the external momenta, and p\ + p 2 + Pz + Pa = 0. Now consider the 
difference 

+ k,p 2 - k,p 3 , Pi ) - T(p 1 ,p 2 ,p 3 ,p i ) = k^(p 1 ,k,p 2 - k,p3,p 4 ) . (37) 

The underlining here refers to the fact that, in the function A M the external 
line with momentum k follows distinct Feynman rules. 

If we follow a path inside the diagram, we can consider the entire ex- 
pression for A M as being built from expressions containing differences. For 
instance, in the propagators: 

i L_ = ^ (-2p-kY 

(p+k) 2 + m 2 p 2 + m 2 [(j> + k) 2 +m 2 ]{p 2 +m 2 )' 1 ' 
or else, in the 3-vertices: 

(p + k) v - p v = k» 6$ . (39) 

We notice that the expressions for A M are all (superficially) ultraviolet 
convergent! Actually, one may set up unambiguous Feynman rules for 
A fJ -(pi,k,p 2 ,—kp3,p4) and observe that this amplitude exactly behaves as 
a 5-point diagram, hence it is (superficially) convergent. 
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For the 3-point and the 2-point diagrams, one can do exactly the same 
thing by differentiating more than once. In practice, what one finds is that 
there is a set of rules containing fundamental irreducible 2-, 3- and 4-point 
vertices, and in addition rules to determine their differences at different val- 
ues of their momenta. The complete procedure thus leads to the following 
situation. 

We start by postulating the so-called 'primary vertex functions'. These 
are not only the irreducible 2-point functions T[ 2 ](p,—p), the irreducible 3- 
point functions ^[3]{pi,P2, — Pi — Pi) and the irreducible 4-point functions 
r[4](pi, • ■ • ,Pi), but also, in addition, the difference functions A™ (p,k, —p—k) 
and A^(pi,k,p 2 — k,—p\ — p 2 ), and finally the functions C/pj, obtained by 
differentiating A™ once more: 

Ara(Pi +tf)P2 - q, -pi -p 2 ) - Af 2] (p 1 ,p 2 ,-pi -p 2 ) = 

(40) 

q v U^(p! - q,g,p 2 , -pi - p 2 ) , 

where one of the other external lines, p\ or p 2 , is underlined. The double 
underlining is here to denote that the two entries are to be treated distinctly 
(because of the factor k^ , the functions are not symmetric under inter- 
change of k and q). 

These primary vertex functions are derived by first considering the differ- 
ences for rj 4 ], A^j and at two different sets of external momenta. These 
expressions are handled as if they were irreducible 5-point diagrams. These 
are expanded in terms of planar diagrams, where all irreducible subraphs of 
4 or less external lines are bundled to form the primary vertex functions. At 
one of the the edges of such a diagram, we then encounter one of the functions 

This way, we arrive at difference equations for the primary vertices, with 
on the r.h.s. again the primary vertices. The primary vertices rp], Aj^ and 
rpi are then obtained by integrating equations such as (37) with respect to 
the external momentum k. The integration constants will be associated with 
the values of the vertices and propagators in the far ultraviolet region, where 
the theory is asymptotically classical. This completes the procedure to obtain 
all amplitudes by iteration. 

Technical implementation of our scheme requires that in all diagrams, an 
unambiguous path can be defined from one external line to another. In QED, 
one may use the paths defined by the electron lines. In a planar theory, one 
may define the paths to run along the edges of a diagram. 

Three remarks are in order: 
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i. The procedure is effectively a renormalization group procedure. The 
functions A M and LP 11 ' play the role of beta functions. 

ii. The procedure is essentially still perturbative, since the planar diagrams 
must still be summed. Our beta functions are free of ultraviolet diver- 
gences, but the summation over planar diagrams may well diverge. 

in. The procedure only works if the integrations do not lead to clashes. This 
implies that it is not to be viewed as a substitute for regularization proce- 
dures such as dimensional regularization. We still need dimensional reg- 
ularization if we want to prove that the method is unambiguous, which, 
of course, it is in the case of planar QCD. 
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